Geometric phases and Wannier functions of Bloch electrons in 1-dimension 
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We present a formal expression for Wannier functions of composite bands of 1-D Bloch electrons 
in terms of parallel-transported Bloch functions and their non-Abelian geometric phases. Spatial 
decay properties of these Wannier functions are studied in the case of simple bands of 1-D model 
insulator and metal. Within first-principles density functional theory, we illustrate the formalism 
through the construction of Wannier functions of polyethylene and polyacetylene. 
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Since the introduction in 1937, Wannier functions 
(WF) 1^ have been conceptually interesting and useful in 
studies of electronic properties of materials. They high- 
light the atomic orbital character of electrons and chem- 
ical bonding 2] and appear naturally in the symmetry 
and analytical properties of energy bands in solids 
They form an excellent basis set for use in the construc- 
tion of first-principles model Hamiltonians in studies of 
phase transitionsj5i,||6|]. There is growing interest in WFs 
as they provide a localized basis in real space for the 
description of electronic states and are useful in the im- 
plementation of computationally efhcient "linear scaling" 
algorithms |3,H- 

WFs are obtained as a unitary transformation of Bloch 
functions 1] . A relationship between the first moment of 
squared WFs (density) and geometric phases j3] of Bloch 
electron was uncovered in the context of adiabatic dy- 
namics of Bloch electron To^ and later in the modern 
theory of polarizationjld^ J^J. WFs are not unique due 
to the freedom in choice of phases accompanying Bloch 
functions, but they are desired to be highly localized for 
practical applications. The second moment (cumulant) 
of squared WFs, a measure of their spatial range, is used 
as an objective function for minimization in the varia- 
tional procedure to construct WFs developed by Marzari 
and Vanderbiltfl^. 

In this paper, we focus on the solution of WF problem 
for systems periodic in 1-D, and comment on its relation 
and extension to higher dimensions. We use Bloch func- 
tions parallel transported along the Bloch vector k in the 
Brillouin zone as a starting point in the construction of 
WFs. As these are not periodic in reciprocal space, we 
present an explicit form for a unitary transformation to 
construct a set that is periodic as a function of k. For the 
1-D case, we show that the WFs obtained as a Fourier 
transform of this set are maximally localized. 

Starting with a subspace Sko of N energy eigenfunc- 
tions at ko, we use parallel transport to obtain Bloch 
functions < 
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.^{x) for k G {ko,ko + ^), a being the lat- 
tice constant and n being an index for electronic states. 



This is accomplished using 



dk 



> obtained routinely 



within the parallel transport gauge: 







(1) 



in density functional theory linear response |14lll5j | (or k.p 
perturbation theory) calculations. We use second-order 
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The electronic part of the macroscopic polarization|llj P 
is readily given by ^^Tr{r). 

Non-diagonal form of the matrix F, in general, in- 
dicates that a starting wavefunction at kg mixes with 
other members of the subspace Stg as k is transported 
to /co + ^ 7 amounting to a unitary transformation of the 
subspace Sk„ ■ Diagonalization of F and corresponding 
transformation of Bloch functions in S^o result in states 
that are decoupled from each other during parallel trans- 
port, effectively reducing the problem of composite bands 
to A'^ simple bands. If M forms a.n N x N matrix whose 
columns are the eigenvectors of F, we obtain an auxiliary 
set of parallel transported wavefunctions Unk ■ 
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where jn^mn — {M^TM)mn- The resulting wave func- 
tions {|Mnfc)} E^re independent of the choice of ko and 
satisfy: (umfeo | exp (z^) |{i„(fe„+277/a)) = exp{ij„)S„in 
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= < X >r 



where 7„ is the n-th eigen value of F and 7„ 

Wavefunctions {\unk)} are not periodic in k with a pe- 
riod ^ if 7„ 7^ 0. A simple U(l) transformation can be 
used to derive {\vnk)} that have the lattice periodicity in 
both reciprocal and real space: 

|u„fc) = exp |-z/c (^^^^ I |u„fc) (4) 
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TABLE L Polarization < a:: > and squared localization length 
P calculated for a single isolated band {n — 1) for different 
potentials: (1) Gaussian potential with c = -~lQ.Q,d — 0.3; 
and (2) asymmetric potential with di — —0.3, d2 ~ —0.2, Ci = 
—5.0, — C2 = —3.0 both with a — 1.0 (all in atomic units). 
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C/„™,fc=exp{-zfc(|^)}M*„ (5) 



In the frame-work of differential geometry|l9| the 
transformation [/„m,fc projects the } in the ray 

space onto {|wn/c) } in the physical space. Each function 
\vnk) is periodic and continuous in fc— space and can be 
Fourier transformed to obtain the WF Wn{x, R): 

k 

- Y^e^^^-^-^'^CnrlUrk) (6) 
k.r 

where are the energy eigenf unctions obtained 

from diagonalization and Cnr — {'Vnk\urk)- Eqn ^ is 
the main formal result of this paper that gives WF in 
the general case of composite bands for systems that are 
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FIG. 1: (A) WF constructed using for the lowest non- 
degenerate band (n=l) in case of ID insulator modeled with: 
Gaussian centro-symmetric potential (solid line) with c — 
— 10, fo — 0.3 and asymmetric potential (dashed line) with 
di = -0.3, d2 = -0.2, ci = -5.0, -C2 = -3.0, both with 
a = 1.0. (B) and (C) shows decay of WF and E(R) re- 
spectively for the Gaussian centro-symmetric potential. The 
power law exponents appear as the coefficients of ln(x) and 
ln(R) in (B) and (C) respectively (all in atomic units). 



periodic in one dimension. Using Eqn and these 
WFs can be shown to have two important properties (in 
addition to being orthonormal) : 

1. J xW*ix,0)Wnix,0)dx=Jn 

2. el,p^„ = J x'w:{x, 0)Wn{x, 0)dx - (7) 

squared localization length i'^yp n independent of phase 
factors accompanying the energy eigenfunctions and can 
be shown to be equal to the one obtained using linear 
response |20l El| . Thus, these WFs are the same as max- 
imally localized Wannier functions|l3|| in 1-dimension. 
A general form of U{1) transformation in Eqn Q is 
exp |— ifc (-5^) — ik.La — ig{k)^, where L is an integer 
and g{k) a periodic function of k. It is readily shown 
that the term k.La in the exponent altersfllj < x >„ by 
La but keeps l^p„ unchanged. On the other hand, the 
term g{k) keeps < x >„ unchanged, but increases 

We propose a generalization of expression © to met- 
als by including the Fermi-Dirac distribution factor 

f{Enk,Ty. 

\W^{x,R)) = ^e*'=(^-^)C„. [f{Erk,T)f- Kk) (8) 

k,r 

These generalized WFs do not form an orthonormal set 
in general, but yield the same form of density matrix for 
metals as that for insulators: 

p'^ix, x') = Wl^x, R)Wl{x\ R). 

nR 

Secondly, they can be shown to have infinite localization 
length in the case of partially filled bands, arising from 
the discontinuity in f{Enk,T) at the Fermi energy '22] . 

We first use above formalism to construct WFs for an 
isolated single band in ID model insulator and metal 
(partially filled band). We consider two generic types 
of model potentials: (1) Gaussian centro-symmetric po- 
tential: V{x) — X^T^i ^j^^C2:p[— (a; — na)^ /d^], which 

was used by He and Vanderbiltji^l in analysis of spa- 
tial decay of WFs and is used to benchmark our calcu- 
lation, and (2) general asymmetric potential: V{x) = 
ci [1 -I- cos {2Tr{x + di)/a}] + C2[l + cos Unix + d2)/a}], 
which we use to test the universality |23| of decay expo- 
nents of WFs. Fig^^A) shows WFs constructed using 
our formalism for the two potentials. 

We calculate the squared localization length using 
three different methods: ^^^„ using Eqn {T)), ^^^^ „ = 
( ^g^" I ^q'^" ) obtained using linear response and i%pn — 

^ T,Bzi^- Re{u^l ,,\u^l i,+^i,))Ak, calculated from the 
parallel transported Bloch functions. As summarized in 
Table H £2 calculated using the three methods for both 
the potentials studied here are identical within numerical 
errors. As the ij^j^ „ is known to be gauge invariant for 
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FIG. 2: Polyethylene: (a) Band structure, (b) Transformed band structure of {unk\H {k)\unk) , (c) WF corresponding to C-H 
bond, (d) WF corresponding to C-C bond. 
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FIG. 3: Polyacetylene: (a) Band structure, (b) Transformed band structure of {unk\H{}i)\unk) , (c) WF corresponding to C-C 
cr— bond, (d) WF corresponding to C-C tt— bond. In (c) and (d) the isosurfaces are prepared using XCrysDen software |T^. 



a given subspace|20j, and correspond to maximum local- 
ization, this numerically verifies maximal localization of 
the WFs constructed using the present formalism for the 
two cases. 

We now compare the spatial decay behavior of the 
present WFs with the maximally localized WFsJ^ which 
exhibit a power law decay in addition to the exponen- 
tial one|2J|: Wn{x) ~ x~'^l'^e}^^ ^ where h is defined from 
the fc— point tt /a-\- ih at which Eq = Ei in the complex 
band structure. Similarly En{R), defined as En{R) — 
{Wn{x,R)\H\Wn{x,{))), decays as: E{R) « x'^/'^e'^^. 
The power law decay exponents —3/4 for Wn{x^R) and 
—3/2 for En{R) are universal numbers[23j. Calculating 
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FIG. 4: Decay of WF for metal with half filled band modeled 
by a centro-symmetric potential at T=OK: dl — 0.25, cl — 
—5.0, ci2 = 0, c2 = with a = 1.0 ( all in atomic units ). 



h from the complex band structure and fitting these gen- 
eral decay forms to the WFs in this work, we find the 
same power law exponents for both the model potentials 
considered. Results for the centro-symmetric potential 
are shown in FigQlB) and (C). 

WF constructed using Eqn. © for a partially filled 
band (T = OK) is shown in Fig01 Form of the spatial 
decay of this WF agrees with that of the density matrix 
(oc a;"'') predicted analytically by Beigi and Arias |2^. 
We observe that the exponent rj is between 0.95 and 1 
for any fraction of band filling at OK. This means the 
localization length of metals being infinity. 

We now illustrate the present method within first- 
principles density functional theory in the general case 
of composite bands by constructing WFs for two poly- 
mer chain molecules: polyethylene (PE) and polyacety- 
lene (PA). A unit cell of a polyethylene (polyacetylene) 
chain consists of 2[CH2] (2[CH]) with twelve (ten) valence 
electrons occupying the lowest six (five) bands. The unit 
cell considered is large enough in the transverse direc- 
tions to minimize the interaction between the periodic 
images of the chain. We use standard density functional 
theory codeji^ to carry out these calculations. We used 
100 Runge-Kutta steps in the parallel transport and 20 
fc— points in the Fourier transform (Eqn JB))) to obtain 
WFs. 

The electronic structure of PE [FigI21^a)] apparently 
groups into two non-crossing sets of bands (the lowest 
two and the rest). Interestingly, unitary matrix M we 
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TABLE II; Polarization < r > and squared localization length 
£wp of polyethylene WFs (all lengths are in A). 



Bond 




<f> 


t-WF 


C-C 


(0.0,0.0,0.0) 


(1.25.0,0.0) 


0.24 


C-H 


(=F0.63,0.53,±0.91) 


(±0.63,-0.53,t0.73) 


0.23 



determined from the geometric phases mixes bands from 
these groups. We use M to transform the Unk to Unk, 
and plot a transformed band structure [Fig|21[b)] with 
expectation values of energy of Unk- We find two dis- 
tinct bands: (1) the lower one with double degeneracy 
corresponding to C-C bonds and (2) the upper one with 
quadruple degeneracy corresponding to C-H bonds. In 
Table ^ we summarize < r > and £^ values for the 
six WFs in two groups corresponding to the two trans- 
formed bands (1) and (2). The average squared local- 
ization length i^p = is 0.23 1^ and 



''WF,n 



agrees well with linear response calculation of value 
0.23 indicating maximal localization. Contour-plots 
of the two WFs are shown in Fig|2Ic,d) 

Similarly, electronic structure of PA (FigO^a)) yields 
the transformed band structure (Figl^Jb)) consisting of 
four isolated bands: two of them (lowest in energy) are 
almost degenerate corresponding to the two C-C cr-bonds 
of slightly different lengths. The one in the middle is 
doubly degenerate and corresponds to the two C-H bonds 
and the one highest in energy corresponds to the C-C 
TT-bond. FigOl^Cjd) shows isosurfaces of WFs describing 
C-C a bond (i^p = 0.221^) and C-C tt bond {l^yp = 
2.95A^) in polyacetylene. An order of magnitude higher 
iyyp of C-C TT bonds than the others makes 7r-conjugated 
systems particularly interesting for transport properties. 

While the logic behind the present formalism is the 
construction of simple continuous and periodic functions 
from parallel transported Bloch functions, the resulting 
WFs in 1-D are maximally localized. Its generalization to 
higher dimensions is natural when the geometric phases 
Ta for paths along different directions in the Brillouin 
zone commute. In most cases, however, the Tq's do not 
commute (corresponding to nontrivial phases for closed 
paths within the Brillouin zone) implying the WFs can 
not be perfectly localized in all the three directions si- 
multaneously. 

In conclusion, we have presented a formal expression 
for WFs of 1-d Bloch electrons and a first-principles 
method for their construction through calculation of ge- 
ometric phase matrices T. It can be readily implemented 
in any DFT code that includes linear response. It natu- 
rally identifies the groups of bands that make up WFs (or 
bonds) through the eigenvectors of the F and their po- 
larization through eigenvalues of F. The method should 
be very useful in studies of 1-D periodic systems in nano- 



science. 
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